function [xArr, qArr] = iterate_nhmm(u, k, a, b, x0, q0, N)
% 多步迭代离散系统(5)
%   x_{n+1} = u sin(x_n)(1 - sin(x_n)) + k(a + b q_n^2) x_n
%   q_{n+1} = q_n + 0.1 x_n
%
%   输入:
%       u, k, a, b  - 系统参数
%       x0, q0      - 初始状态
%       N           - 迭代次数
%   输出:
%       xArr, qArr  - 分别存放每一步的 x_n 和 q_n

    xArr = zeros(1, N);
    qArr = zeros(1, N);

    xArr(1) = x0;
    qArr(1) = q0;

    for n = 1:(N-1)
        x_n = xArr(n);
        q_n = qArr(n);
        [x_next, q_next] = nhmm(x_n, q_n, u, k, a, b); % 调用您的 nhmm.m
        xArr(n+1) = x_next;
        qArr(n+1) = q_next;
    end
end
